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The selective optical detection of individual metallic nanoparticles (NPs) with high spatial and 
temporal resolution is a challenging endeavour, yet is key to the understanding of their optical re- 
sponse and their exploitation in applications from miniaturised optoelectronics and sensors to med- 
ical diagnostics and therapeutics. However, only few reports on ultrafast pump-probe spectroscopy 
on single small metallic NPs are available to date. Here, we demonstrate a novel phase-sensitive four- 
wave mixing (FWM) microscopy in heterodyne detection to resolve for the first time the ultrafast 
changes of real and imaginary part of the dielectric function of single small (< 40 nm) spherical gold 
NPs. The results are quantitatively described via the transient electron temperature and density in 
gold considering both intraband and interband transitions at the surface plasmon resonance. This 
novel microscopy technique enables background-free detection of the complex susceptibility change 
even in highly scattering environments and can be readily applied to any metal nanostructure. 



Metallic nanoparticles exhibit morphology-dependent 
electromagnetic resonances also called surface plasmon 
resonances (SPR) which couple to propagating light. 
These resonances originate from a coherent oscillation 
of electrons in the metal where the restoring force is due 
to the electric field created by the corresponding charge 
displacement. Since this electric field depends on the 
morphology, the SPR frequency and linewidth is a func- 
tion of the shape and size of the NP and its dielectric 
environment—"—. The resulting local optical resonances 
can be exploited to image metallic NPs with high spa- 
tial resolution and to probe nanoscale regions in the NP 
vicinity through the local electric field of the resonance. 
Especially gold NPs are ideal optical labels for biological 
applications owing to their bio-compatibility and photo- 
stability, and much effort has been devoted recently to 
develop techniques capable of detecting gold NPs in cells 
and tissues with high contrast and sensitivity to the sin- 
gle particle level—"—. 

Ultrafast optical spectroscopy of metallic NPs is an 
intriguing area of research investigating the correlated 
electronic and vibrational dynamics in nanosized met- 
als. Due to limited detection sensitivity, most of the 
experiments performed so far— lii used large ensembles 
of metallic NPs, with the associated difficulty of inho- 
mogcncous NP size, shape and environment fiuctuations. 
Only few reports on ultrafast pump-probe spectroscopy 
on single small metallic NPs—"— are available to date. 
Importantly, these experiments lacked the ability to ex- 
tract the change in the NP dielectric function as a com- 
plex quantity separating its real and imaginary part—. 
This however is the key physical quantity of interest, be- 
ing a function of the electron and lattice temperatures 
transiently changing during the thermalization dynam- 
icsiS. 

To overcome these limitations we have developed a 
phase-sensitive transient four-wave mixing microscopy 
technique as sketched in Fig.[T] Laser pulses of ~ 150 fs 
Fourier-limited duration with a tuneable center wave- 



length from 540 nm to 590 nm and i^l ~ 76.1MHz repe- 
tition rate are split into three beams. One beam acts 
as pump and excites the NP with an intensity which 
is temporally modulated by an acousto-optic modula- 
tor (AOM) driven with an amplitude modulation of fre- 
quency t'm = IMHz. This implementation replaces the 
two beam interference used in our previous work—, pro- 
viding a stable phase of the modulation. The change in 
the NP optical properties induced by this excitation are 
probed by a second pulse at an adjustable delay time 
T after the pump pulse. A FWM signal proportional 
to E1EIE2 with El, E2 electric fields of the pump and 
probe, respectively, is collected in transmission and de- 
tected interferometrically using the third beam acting as 
reference. We used a heterodyne scheme to discriminate 
the FWM field from pump and probe pulses and to de- 
tect amplitude and phase of the field. By up-shifting the 
probe optical frequency via a second AOM (driven with 
a constant amplitude at a frequency of V2 — 78.8 MHz) 
we detected the interference of the FWM with the ref- 
erence at the frequency 1^2 — ^^m ^ t'L = 1.7 MHz by a 
dual-channel lock-in amplifier. We simultaneously detect 
the interference of the transmitted probe with the refer- 
ence at V2 — v-L = 2.7 MHz using a second dual-channel 
lock-in, and deduced the amplitude {A) and phase (P) 
components of the FWM field relative to the transmitted 
probe (see sketch in Fig.[T] and supplementary material 
(SM)). We have previously demonstrated background- 
free FWM imaging in a highly scattering and fiuoresc- 
ing environment, and sensitivity to the single NP level 
with a spatial resolution significantly surpassing one pho- 
ton diffraction—. An example of FWM imaging of gold 
NPs of 20 nm diameter resonantly excited and probed at 
550 nm is shown in the inset of Fig. [1] together with an 
histogram indicating a monomodal size distribution well 
separated from the noise level (here ~ 1/100 of the NP 
signal), confirming that each spot corresponds to a single 
NP. 

To explain how the detection of A and P enables us to 
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FIG. 1. Sketch of the phase-resolved transient four-wave mixing micro-spectroscopy technique. A pump pulse train excites 
the gold NP SPR with a temporally modulated intensity. The resulting change in the NP optical properties is probed by a 
second pulse train frequency shifted by an acousto-optic modulator (A0M2) and the FWM field is collected in transmission and 
detected interferometrically at the corresponding mixed frequency using an unshifted reference. BS: 50:50 beam splitter, BPD: 
balanced photodiodes. Top inset: FWM imaging of 20 nm gold NPs resonantly excited and probed at 550 nm and corresponding 
histogram with a Gaussian fit showing a monomodal size distribution (the deduced relative size distribution is 30%). Right: 
Sketch of the phase-resolved detection using two dual-channel lock-in amplifiers to measure the X and Y components of the 
probe field E2 and FWM field -Efwm- 



distinguish the transient changes of the real and imag- 
inary part of the dielectric function of a single NP, we 
write the induced polarization field at the particle as 
P = eofdQ;i?2 where eo is the vacuum permittivity, ed is 
the dielectric constant of the medium surrounding the NP 
and a is the particle polarizability. For a particle radius 
R much smaller than the wavelength of light (Rayleigh 
limit), the polarizability is given by a = 47r_R'^ l+2t^ with 
the particle dielectric constant e. We denote ed 
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as an effective dielectric constant, and the FWM field 
in terms of the nonlinear polarization from the pump- 
induced change Ae^, ie -Efwm oc AeE2. The complex 
dielectric constant is then expressed in terms of its ampli- 
tude and phase Ae'= |Ae] exp(i(/?) — Aer -I- iAej, as well 
as the probe field E2 = |i?2| exp(i^), such that i?FWM 
|A?||i?2| exp(i((/) -I- (/?)). Hence (p is the phase difference 
between the FWM field and the probe field at the particle 
position. The transmitted probe in the far field acquires 
a phase shift compared to a spherical wave of the FWM 
field, known as Gouy phase shift for a Gaussian beam. 
The Gouy phase shift is — 7r/2 from the focus to infinity, 
hence in the far-field the FWM field is phase-shifted rela- 
tive to the transmitted probe by ip + Tr/2. The amplitude 
modulation of the transmitted probe is then (see sketch 
in Fig.HI) A cx |A?1 cos((^ -I- tt/2) = -|Ae| sm{ip) ^ -AeJ 
and the phase modulation is P oc |Ae|sin((p + 7r/2) — 
|Ae]cos((/p) = Aer. Note that the particle absorption 
cross-section is proportional to the imaginary part of a 
as can be deduced from Mie theory in the small particle 
limit, hence ^ is a measure of the change in absorption 



cross-section. 

The measured A and P versus pump-probe delay for 
single gold NPs of various sizes resonantly excited and 
probed at 550 nm are shown in Fig. [51 Both exhibit an 
initial dynamics in the picosecond time scale followed by 
a decay in the hundreds of picoseconds. Qualitatively, 
one can understand these dynamics as being related to 
the pump-induced increase in the electron temperature 
which provokes a broadening and a shift of the SPR 
probed by £2- The picosecond decay is related to the 
thermalization of the hot electron gas with the cooler 
lattice by electron-phonon coupling, while the thermal- 
ization of the heated NP with the cooler surrounding is 
providing the subsequent decay. To achieve a quantita- 
tive description of the measured dynamics we modeled e 
using the gold NP dielectric function 
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where ujp = Uoe'^ / eQirio is the plasma frequency (ue, e and 
mo being the conduction electron density, charge and ef- 
fective mass respectively), F is the inverse of the electron 
relaxation time, and is the contribution due to bound 
electrons associated with interband transitions from the 
d-bands to the conduction band. At room temperature in 
the absence of optical excitation F can be approximated 
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where 70 is the bulk damping rate, wp = 1-4 x lO^m/s 
is the Fermi velocity, and g parametrizes surface damp- 
ing effects. For bulk gold we used fiijj-p = 9eV and 
^7o = 70meV, well reproducing the Drude part of 
the bulk dielectric constant measured by Johnson and 
Christy These values also reproduce experimental ab- 
sorption spectra of single gold nanoparticles in the liter- 
ature™. The parameter g varies in the literature^''' from 
0.2 to 2.2, and even between individual NPs within the 
same experiment. This variance is likely to be related to 
the atomistic NP surface structure/faceting— hence we 
use g as free parameter falling in the range of previously 
reported values. In the presence of optical excitation 
we then model the change in the electron temperature 
To and lattice temperature T\ using a two-temperature 
model^^ 

Cc{T,)^ = -Gcp(re - TO + Pit), (3) 



a^ = Gep(Te-ri)-g, (4) 

where Co is the electronic heat capacity, well described 
by Co = for^s To < 3000 K, Ci is the lattice heat ca- 
pacity, and Cop is the electron-phonon coupling constant. 
Co and C\ are taken as the bulk value scaled to the NP 
volume (see SM) . Q describes the heat loss from particle 
to the surrounding and is calculated using the diffusion 
equation for the temperature of the surrounding medium 
(see SM). P{t) describes the excitation process which has 
two components, P{ and Pb, related to free and bound 
electrons, respectively. The bound electron contribution 
has to be considered since in spherical gold NPs the onset 
of interband transitions from d-bands to the conduction 
band at the Fermi surface is energetically close to the 
SPRj^ Pf is modeled as a rapid (~ 10 fs, quasi instanta- 
neous for the 150 fs pulses used) transfer of energy from 
the resonantly excited SPR to single electron excitations 
of the conduction band, which thermalize with the free- 
electron gas creating a hot Fermi-Dirac distribution on a 
characteristic time scale rt— li^. Similarly Pb originates 
from the transfer of energy from the SPR to interband 
excitations and consequent heating of the free electron 
gas via Auger electron-hole recombination with a char- 
acteristic time Tb (see SM) . Noticeably, most pump-probe 
experiments in previous literature use a pump photon en- 
ergy below the onset of interband transitions, and thus 
neglect the latter effect— lii. The time integral of P(t) 
is equal to the absorbed energy calculated from the NP 
absorption cross-section and the pump pulse fluence in 
the experiment. 

The non-equilibrium electron and lattice temperatures 
modify the electron-electron and electron-phonon scat- 
tering processes, which can be described^- by F = Fq + 
7op + 7co with 

7cp(?i) = bcph^uj'^iTi - To), (5) 
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FIG. 2. Transient changes of the real (Aea) and imaginary 
(Aei) part of the dielectric function of a single gold NP res- 
onantly excited and probed at 550 nm, for different NP di- 
ameters as indicated. Pump (probe) fluence is 0.60 J/m^ 
(O.llJ/m^) for the 15 nm and 20 nm NPs and 0.65 J/m^ 
(0.05 J/m^) for the 38 nm particle. Acquisition time per point 
is 200 ms. Dashed hnes are corresponding calculations (see 
text), with the green lines neglecting the heating by interband 
absorption. Curves are vertically displaced for clarity. The 
top sketch illustrates the SPR excitation by the pump (Ei) 
and the subsequent heating and cooling dynamics monitored 
by the probe. The inset shows coherent phonon oscillations 
for the 38 nm NP. 

and 

7oo(To) -6oo7iV(T2-T2). (6) 

using the equilibrium temperature Tq. The electron and 
lattice temperatures also affect e^. To estimate this effect 
we calculated the imaginary part €2 of due to inter- 
band transitions near the X and L points, starting from 
the description in Refs.'20' and[2ll. The temperature de- 
pendence of e** results from the Fermi distribution and a 
spectral broadening proportional to F (see SM). The real 
part of e** is calculated from €2 using Kramers-Kronig. 

The transient change in e modifies ^ and is the source 
of the measured AeJ and Aer via A = — ^Ae} and 
P = ^AeR (see Fig. (2) with a constant ^— . We are able 
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FIG. 3. Transients of the electron temperature (solid lines) 
and lattice temperature (dashed lines) for different NP sizes 
as indicated, calculated using Eq.|31l3] for a pump fluence of 
0.60 J/m^. The dotted curve is the electron temperature of a 
20 nm NP when the excitation term involving bound electrons 
is neglected. The inset shows the electron-phonon coupling 
parameter used in the calculations. The dashed line in the 
inset indicates the volumetric scaling. 



to consistently model all measured dynamics for particle 
sizes up to 40 nm diameter and for different pump exci- 
tation intensities in the 0.16 — 1.3 J/m^ range with only 
few particle-dependent (g, Gcp) and particle independent 
(6oo, &cpj Tf, Tb) parameters. The gold band-structure 
parameters to calculate are taken from literature— i^i. 
The proportionality term ^ scales with and is indepen- 
dent of the pump intensity as expected from the particle 
polarizability a. Larger particles have a polarizability 
which deviates significantly from the Rayleigh limit used 
here, and we found that FWM dynamics measured on 
100 nm diameter particles can no longer be described by 
our model (see SM) . We calibrated the size of the 38 nm 
NP investigated in Fig. [5] using the period of the coherent 
phonon oscillations observed in AeR (see inset in Fig. [2]). 
These are due to the modulation of the plasma frequency 
by the breathing vibrational mode of the nanoparticle 
with an oscillation period proportional to the particle 
radius^iHii^(see SM). On NPs smaller than 20 nm the 
oscillation period was not resolved, possibly because its 
time scale becomes comparable to the lattice expansion 
which is no longer impulsively exciting the vibrational 
mode hence its amplitude falls below the measurement 
noise. In this case the particle size was derived from the 
strength of the measured FWM field relative to that of 
the 38 nm particle, consistently with ^ oc i?"^ . 

We use Tf — 500 fs independent on particle size in 
the 15-40 nm diameter range from Ref.lm bcc = 0.89 x 
lO^s-^eV-^K-^ and 6op = 0.5 x lO" s-^eV^K-i con- 
sistent with Ref.flQl, and g in the range 1.8-1.9 consistent 
with Ref.i. 

The fast rise of AeJ and AeR observed in our data re- 
veals a fast d-state Auger recombination time of Tb = 
70 fs ^ Tf. While this time has not been measured pre- 



viously in Au, it is within the energy-dependent range of 
250-20 fs measured in Cu— . To clarify the significance 
of interband transitions in the SPR resonant excitation, 
we show in Fig. [2] also simulations (green lines) neglecting 
the heating Pb in Eq.[3]due to the interband transitions. 
The effect of interband transitions in the excitation is 
important to account not only for the magnitude of the 
measured change, but also its dynamics, particularly for 
the observed fast rise time. We could simulate the mea- 
sured rise time without the interband contribution only 
when using Tf < 100 fs, much faster than the well charac- 
terized 500 fs electron thermalization dynamics-^ in bulk 
gold and gold NPs larger than 10 nm. The faster Auger 
recombination of the d-holes thus gives rise to a quali- 
tatively different initial dynamics, and can be attributed 
to a higher density of states available in the scattering 
involving the lower-lying d-states. 

It is instructive to examine the temperature transients 
of To and T\ plotted in Fig.[3l Noticeably the electron 
temperature reaches a maximum value which is almost 
independent of NP size. This is a consequence of the vol- 
umetric dependence of both the electron heat capacity 
and the absorption cross-section, balancing each other. 
The decay of To on the picosecond time scale coincides 
with the rise of Ti due to electron-phonon coupling, and 
the electron temperature after thermalization is much 
lower due to the significantly larger lattice heat capacity 
(Ci >• Co). We find that Gep also has a nearly volu- 
metric dependence, and in turn the electron-phonon dy- 
namics exhibits only a weak size dependence in the size 
range 15-40 nm investigated here, consistent with reports 
in the literature We also observe that the cooling with 
the environment is size dependent, with smaller particles 
cooling faster^ ^. We point out that the average temper- 
ature increase at the particle surface is only 1.3 K, 1.9 K, 
and 4.7 K for 15 nm, 20 nm and 38 nm diameters respec- 
tively and can be lowered to 1.3 K for 38 nm diameter 
at the lowest (0.16 J/m^) pump excitation used in the 
experiment. Importantly, such low average photother- 
mal heating indicates that our technique is compatible 
with live cell applications. Finally, the effect of the inter- 
band excitation Pb is shown by plotting Te in its absence 
(dotted curve in Fig. [3]), which shows a delayed rise and 
reaches only about half of the maximum electron tem- 
perature. 

To visualize the origin of the FWM signal in terms of 
broadening and shift of the SPR we plot in Fig. [4] the 
spectra of ei and cr at different delay times for an exci- 
tation at 550 nm. The negative change of ei probed at 
550 nm (see dotted red vertical line in Fig. 3]) is mainly a 
result of the broadening of the SPR and in turn quench- 
ing of the resonance. Conversely, the negative change of 

at 550 nm is a combination of the SPR broadening and 
its frequency shift toward lower energies. Measurements 
and simulations of Ae} and AeR excited and probed at 
540 nm {Tiuj = 2.3 eV) and 575 nm {huj = 2.16eV) eluci- 
date the cases of AeJ being mostly sensitive to the SPR 
broadening and Acr to the SPR shift, and vice versa. 
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FIG. 4. Left: Spectrally-resolved real and imaginary part of the dielectric constant e simulated for the 38 nm NP at various 
delay times, for an excitation at 550 nm. Dashed lines are the differential spectra at r = 0.5 ps. Vertical dotted lines indicate 
the excitation and detection photon energies used in the experiment on the right. In the inset, the phase of e is shown at 
r = 0.5 ps (symbols: experiment, line: simulation). Right: Measurement and simulations of the transient changes of the real 
and imaginary part of the dielectric function at different photon energies, as indicated. The pump fluence was 0.48 J/m^ at 
540 nm and 0.6 J/m^ for all other wavelengths. Curves are vertically displaced for clarity. 



Specifically we see that at 2.3 eV, AeR is a sensitive 
probe of the modulation of the SPR by the breathing 
vibrational mode, while Aef is essentially unaffected by 
it. We also performed measurements and simulations at 
590 nm {Tiuj = 2.10 eV) where Aej changes sign. Notice- 
ably, the initial dynamics measured at different wave- 
lengths is consistently reproduced by our model of the 
fast Auger recombination of the d-holes. Finally, to il- 
lustrate that the novel imaging technique presented here 
can be used to sense SPR shifts relative to the probe 
wavelength in an intrinsic ratiometric way independent 
of the signal strength, we show in the inset of Fig. 2] the 
phase If of Ae — \Ae\exp(iip) from the ratio between 
its imaginary and real part at r = 0.5 ps. ip changes by 
nearly 180 degrees over only 50 nm wavelength as a direct 
manifestation of the transition from Aer = 0, Aei < 
{ip = 37r/2, at 540 nm) to AeJ = 0, Ae^ < = tt, at 
575 nm), to AeJ > 0, Aej; < {n/2 < < tt, at 590 nm). 

Our technique is applicable to any NP shape to re- 
veal new physics insights in the electron and lattice dy- 
namics of novel metallic nanostructures. Importantly 
it is compatible with living cell applications and offers 
background-free time-resolved detection of the full com- 
plex change of the dielectric constant. The intrinsic ra- 
tiometric information in the signal phase enables sens- 
ing applications such as monitoring nanoscale distance 
changes with plasmon rulers— in highly scattering envi- 
ronments such as cells and tissues. 

I. METHODS 

Commercially available gold NPs in aqueous suspen- 
sion (BB International) were drop cast on a glass cover 



slip at sufficiently low density so that regions with well 
separated individual NPs were formed, and subsequently 
embedded in a mounting medium (Cargille Meltmount) 
index-matched to the oil-immersion microscope objec- 
tives used for high-resolution imaging (see also Ref.Q). 
The sample is scanned using an xyz piezoelectric stage 
(Physik Instrumente Nanocube). 

Optical pulses were provided by the intra-cavity 
frequency-doubled optical parametric oscillator (APE- 
PP2) pumped by a TiiSapphire oscillator (Coherent 
Mira). Pulses were linearly polarized vertically (V) or 
horizontally (H) in the laboratory system using A/4 and 
A/2 waveplates, and a VHH polarization configuration 
for El, E2 and the reference field was used to suppress 
background from detector nonlinearities— . 
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I. SETUP DESCRIPTION 

Optical pulses are provided by the intra-cavity frequency doubled signal of an optical 
parametric oscillator (APE-PP2) pumped by a femtosecond Ti: Sapphire oscillator (Coherent 
Mira) oi ui, — 76.1MHz repetition rate. We used pulses with a tuneable center wavelength 
from 540 nm to 590 nm and a Fourier-limited duration ranging from 150 fs to 180 fs (intensity 
full- width at half-maximum) depending on the wavelength (corresponding to 2.5 nm spectral 
width at 550 nm). The pulse train at the output of the optical parametric oscillator is split 
into three beams having all the same center wavelength (degenerate FWM scheme). The 
delay times between the pulses are controlled with 10 fs accuracy by hnear stages. Pump and 
probe pulses of fields Ei and E2 respectively are recombined into the same spatial mode and 
focused onto the sample by an oil-immersion microscope objective (MO) of 1.25 numerical 
aperture (NA) to achieve high spatial resolution. The sample is positioned with respect to 
the focal volume of the objective by a xyz piezoelectric stage with 10 nm resolution (Physik 
Instrumente Nanocube). The FWM signal is collected in transmission by a second MO 
opposite to the sample, recombined in a beamsplitter with a reference pulse of adjustable 
delay, and the resulting interference is detected by balanced silicon photodiodes and a lock-in 
amplifier (SRS SR844). Discrimination of the FWM signal is achieved using a heterodyne 
scheme. In this scheme two acousto-optic modulators (AOMs) are used to upshift the optical 
frequencies of £^1,2 by the radiofrequencies ui = 80 MHz and 1^2 = 78.8 MHz, respectively. 
The intensity of the first pulse is modulated with close to unity depth by applying a z/m = 
IMHz square wave amplitude modulation to the first AOM. The interference of the FWM 
signal with the unshifted reference field gives rise to a beat note at the radiofrequency 

— and its sidebands separated by the repetition rate i^l of the pulse train. Using 
an electronically synchronized reference created by mixing the AOM drive frequencies with 
the laser repetition rate, we detect the interference of the FWM with the reference by a 
dual-channel lock-in amplifier at V2 — — = 1.7 MHz. Simultaneously, we detect the 
interference of the E2 field with the reference using a second lock-in amplifier at V2~^h- We 
employed a cross-polarization scheme to reject background originating from nonlinearities of 
the photodiodes, in which the polarization state of Ei was adjusted by a combination of A/4 
and A/2 waveplates to be cross-hnearly polarized to E2 and the reference beam. Polarizers 
in front of the photodiodes transmit the FWM, E2, and reference fields, but block the pump 
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field El. 



II. AMPLITUDE AND PHASE PARTS OF THE FWM FIELD 

Two dual-channel lock-in amplifiers are used to measure the X and Y components of E2 
(probe field) and -Efwm (FWM field), relative to the reference field, as shown in the paper 
in Fig. 1. They are electronically referenced to the frequencies z/2 — i^l and 1^2 — ^'m ~" '^l? 
respectively, in a phase-stable manner. We define the amplitude part A as the projection of 
-EpwM onto E2 and the phase part P as its orthogonal component: 



_ X2XFWM + ^2^FWM 



E,r- « 



We calibrate the relative reference phase of the two lock-in amplifiers using a known 
amplitude and phase modulation of the probe. For this purpose we measure an electrically 
pumped quantum-dot optical amplifier, which shows pump-induced changes of the gain 
and refractive index ^ resulting in large changes of amplitude (> 10%) and phase (> 10°) 
of the transmitted probe field. These changes can be measured with the FWM experi- 
ment, but also directly by measuring X2 and I2 with and without the pump using a slow 
(5 Hz) modulation. Both measurements have to yield the same result, since the response 
of the optical amplifier is the same at 1 MHz and 5 Hz. This requirement determines the 
phase offset 0o between the two lock- in amplifiers in the FWM experiment, which is taken 
into account by transforming the values Xg^"'^'', 1^2*^^^^ directly measured by the lock- in into 
X2 = ^2^^^ cos(0o) — ¥2^^^ sin(0o), Y2 = ¥2'^^ cos(0o) + -^2^"* sin(0o)- The phase offset repro- 
ducibility during the course of the experiments was found to be a few degrees (corresponding 
to a few nanoseconds delay), as expected from the RF electronics used to create the lock- in 
reference signals. 
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III. MODEL OF THE RESPONSE 
A. Excitation 

The response of the GNP polarizabihty to an optical excitation is dominated by the dy- 
namics of the electron distribution. The optical excitation creates an internal optical electric 
field in the material according to the morphology and dielectric material response. This 
field is absorbed creating non-thermal electron distributions, which thermalize by electron- 
electron scattering with characteristic thermalization times, establishing a hot Fermi-Dirac 
distribution with broadened Fermi-edge^i^. 

We can distinguish two different electronic excitations creating a heating power P for the 
electron distribution, related to free (Pf) and bound (Pb) electrons. The latter is important 
in our experiments since the onset of interband transitions from d-bands to the conduction 
band at the Fermi surface is close to the SPR in spherical gold NPs-. We emphasize that 
the dynamics from a single particle excitation to a thermalized electron distribution by 
cascaded electron-electron scattering processes is non-trivial, and depends on the electronic 
band-structure and the initial single particle excitation. It is thus expected that different 
initial single particle excitations are described by different effective thermalization times r. 

The total energy absorbed by the NP is given by /(Jabs, where / is the pump pulse fluence 
and (Tabs is the NP absorption cross section. It is distributed among the different single 
particle excitations (SPE) related to the free electron gas (SPE = f) and the interband 
transitions involving the X and L point (SPE = X, L) respectively. The fraction absorbed 
by the corresponding SPE is given by the relative contribution /spe to the imaginary part 
of the bulk dielectric constant of the NP e2(wexc) at the excitation frequency cjexc, 



where €2 is the contribution in 62 related to the SPE. 

The SPR lifetimes are typically in the ~ 10 fs range and can be neglected compared to 
the pulse duration used in our experiments. Assuming a Gaussian intensity profile of the 
pump pulse 



and a subsequent thermalization described as a Markovian process with time constant r, we 




(2) 




(3) 
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find the heating powers 



r> -^'^abs/sPE 

-rsPEirj = — 

^TSPE 



1 + eri 



- 4trc 



2rsPE0- 

where a is given by the pulse duration in the experiment measured for every wavelength (eg 
a = 0.11 ps for an intensity full- width at half- maximum of 0.18 ps). 

Pf originates from the intraband absorption of the free electron gas and has been described 
in the literature extensively, having a thermalization time Tf ^ 500 fs. 

Pb = Px + Pl instead originates from the absorption of the SPR by interband transitions 
and consequent Auger electron-hole recombination and heating of the free electron gas, 
which we model with a separate thermalization time Tb = tl = tx- 

To calculate the transient hole densities which influence the dielectric response, we express 

P,it) = {nt{t)+<\t))^^ (5) 

in terms of the densities of photoexcited holes n^^, near the X and L points, respectively, 
which lead to Auger recombination heating of the electron distribution. Here V is the particle 
volume. The densities evolve according to 

at Vvror^ V / 
which represents the generation by the Gaussian pulse excitation and the decay due to Auger 
recombination. The photogenerated densities nx,L are given by 

nx,h = -TTT . (7 

At the resonant excitation wavelength of 550 nm in our experiment we find (see section lTlI Cp 
fx = 0.59 and /l = 0.16, i.e. 75% of the excitation is absorbed by interband processes. This 
corresponds for a 38 nm diameter NP with a pump fluence of / = 0.32 J/m^ to photogen- 
erated hole densities of rix = 3 x lO^^rie, riL = 8 x 10~^?2e, with the free electron density 
rig = 5.9 X lO^^m"'^. For the other excitation wavelengths used we find fx = 0.57 and 
/l = 0.24 at 540nm, fx = 0.59 and /l = 0.09 at 575 nm, fx = 0.57 and /l = 0.06 at 
590 nm. 



B. Two-temperature model 



To describe the energy exchange between the optical excitation, electrons in the NP 
and lattice we use a two-temperature model including heat diffusion with the surrounding 
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medium-1^ 



Cj^ = -G^^iT, - Ti) + P,it) + P^it), (8) 



or, ^ 1 d 



dr 



(9) 

r=R 



(10) 



dt (9r 

where = 'jT^ is the temperature dependent electronic heat capacity, Ci is the lattice 
heat capacity, G^p is the electron-phonon coupling constant. We use 'y/V = 67.6 Jm~'^K~^ 
and Ci/V = 2.49 x 10^Jm~^K^^ according to literature-i^, and deduce the particle values 
by multiplying these volumetric constants with the particle volume V = ^nR^. Ts(r, t) 
is the temperature of the surrounding which obeys the radial heat diffusion Eqn.([TO|). 
Ks = 0.07Wm~^K~^ is the thermal conductivity and = 0.1 x 10~^m^s~^ is the ther- 
mal diffusivity of the surrounding. Those values are obtained from the best fit to the FWM 
dynamics. The thermal conductivity and diffusivity of the mounting medium surrounding 
the nanoparticles in our experiment (Cargille Meltmount) are not known, but the fitted 
values are similar to the known values of PVC and epoxy resins. The initial and bound- 
ary conditions are as follows: Electrons, lattice, and surrounding medium are initially at 
thermal equilibrium, ie To(0) = ri(0) = Ts(r, 0) = Tq = 300 K. The surrounding medium 
in contact with the particle surface and the metal lattice are always at thermal equilib- 
rium, Ti(t) = Ts{R,t). The surrounding medium far away from the particle is at ambient 
Ts(oo,t) = 300 K temperature. 

The system of differential equations was solved in MATLAB using the function pdepe, 
which solves initial-boundary value problems for parabolic-elliptic partial differential equa- 
tions in one dimension. 



C. Contribution of interband transitions to the bulk dielectric constant 



To calculate €3 we use the approach developed in Refs. 
the joint density of states (Eqn.(9) in Ref.|7|) 



.^andfl 



We write €2 in terms of 
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where Ax is proportional to the square of the transition dipole matrix element Px of the 
interband transition from the d-band to the Fermi surface near the X point of the band- 
structure. Similarly Aj^+ and are proportional to the square of the transition dipole 
matrix elements near the L point, involving the 5+6 and 4 d-bands respectively (see band 
structure calculated in Ref.jS,). Jx, ■A.+ .g' "^lJ corresponding transition density of 

states which have a temperature dependence given by 



HE) 



{UE,%)-U{E,T,))dE 



(12) 



with i = X, L4 



/p and /d are the occupation probabilities of the p-band and d-band respectively 



l + exp( 



/d(i^,Te) = l 



i -h exp ^^^^ 



and 



\ 



m 



hu- E - hu^. - huo + '-^{E + hoj^.) 



(13) 
(14) 

(15) 
(16) 



where hu is the photon energy, Tiuq is ^Wg^g or huj^+, and correspondingly m^^ refers to 
either the 5+6 or the 4 d-bands. The Fermi level as a function of electron temperature 
Ep{Te) is calculated from the relationship 



DOS{E)dE 



(17) 



-l + exp(^) 

consideririg the total electron density A^^e and using the calculated density of states DOS(-E) 
from Ref.y. To calculate the change on the Fermi Level, we can limit the integration range 
from Ep = Ef{0) - 10 eV to E^ = Ef{0) + 20 eV, 



Sf(0) 



DOS{E)dE 



t DOSiE)dE 
1 + exp 



(18) 
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The limits of integrations in Eq. 
— 20/cbT, and 



are calculated as in Ref. 



.y and B 



and are E^:, 



X 



hoOf^ + {hul — huj, 



X 



m 



m 



X 



m 



X 



if hu < UlOq + ^CJg, or 



E, 



X 



X 



m 



X 

d± 



m5i +m^j_ 



(19) 



(20) 



if hu > UuJq + huj'^Q. Near the L point they are 



if huj < huQ, or 



if hu > huo, and 



E' 



Et 



-hoo^- + {hoo — Tiuq)- 



'd\\ 



m. 



-huj^- + {hco — hojQ 



m 



+ mi 



(21) 



(22) 



E, 



L 

min 



X 




(23) 



"Pll ■ d\\ 

The effective mass parameters and band gap parameters are taken from the Au band 
structure near the X and L points-"^ and are = 0.12 mo, = 0.91 mo, = 0.22 mo, 



m 



X 
dl. 



r + 

^5+6 



0.75 mo, mp|| = 0.25 mo, m^^ = 0.23 mo, rn^jj = 0.63 mo, m 



r + 

^5+6 



0.57 mo, m^l = 0.49mo, hu^ = 1.77 eV, hu^ = 1.47eV, hu^^ = 0.72eV, /iwf+g = l-45eV, 



0.7 mo, m^ii 

^5+6 ~ ^ 



hu^+ = 2.6 eV, and the lattice parameter is a = 4.08 x 10~^°m. 

To fit the measured €2 at 300 K by Johnson and Christy— we use the constants Ax, , 
as fitting parameters similar to Ref.Q, and obtain \Px/ Pj^+^J'^ = 0.26, |Px/-Pl+P ~ 
0.67. One also needs to account for a broadening of eg in order to fit the low-energy tail 
experimentally observed. Instead of introducing an effective temperature as fit parameter-, 
we consider the Drude broadening by convoluting the expression calculated according to 
Eq. [11] with the function 



B{huj) 



1.76 
26f 



cosh 



-2 



'l.76huj\ 



br 



- 



(24) 



which has an integral normalized to 1 and a full-width at half-maximum bT proportional to 
the Drude broadening. We use 



br = 5n 



7o 



+ (6ee(Te' - T^) + b,p{T, - To)) {hu. 



(25) 
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which fits well the data of Johnson and Christy— at 300 K and accounts for additional 
broadenings due to electron-electron and electron-phonon scattering processes in the pres- 
ence of a non-equilibrium electron and lattice temperature on the same ground as for the 
Drude part. 



D. Coherent phonon oscillations 

The visibility of coherent phonon oscillations in Ae^ is described by the change of the 
plasma frequency due to the changing volume of the NP during the breathing oscillation, 
driven by the non-adiabatic thermal expansion of the NP during the thermalization with 
the lattice in the first picoseconds: 



p 



«(Ti-ro) (l + Ccos(^ + 



(26) 



with the bulk volumetric expansion coefficient a = 4.5x 10~^ K^^. To fit the oscillations mea- 
sured on the 38 nm NP shown in the paper (inset of Fig. 2) we used C = 0.75, To = 12.8 ps, 
Td = 80 ps, and = —3.33. The oscillation period is the period of the fundamental breath- 
ing mode in a spherically symmetric NP— , from which we deduce a particle diameter of 



38 nm following the values in Ref.ll2l The quantitative agreement of the oscillation ampli- 
tude using a factor of C = 0.75 is supporting this interpretation, since for an instantaneous 
lattice heating we would expect C = 1, and for a realistic heating time of a few picoseconds 
this value is expected to be somewhat reduced. 



IV. FWM DYNAMICS FOR LARGE NANOPARTICLES 

The SPR of a GNP in a dielectric surrounding can be described by the dipole approxi- 
mation for a GNP diameters of less than about 40 nm. For larger sizes retardation effects 
become significant, resulting in a redshift of the extinction peak with increasing size. To 
verify the suitability of our model to describe large NP sizes, we have investigated the dy- 
namics of Ae for a 125 nm GNP excited and probed at 550 nm (solid lines in Fig.[T]). The 
size of the NP has been estimated by the oscillation period measured in the phase part of 
the signal P (see previous section). A major difference with respect to the smaller NP sizes 
investigated is that P has a positive sign. Within our model in the dipole approximation, 
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FIG. I. Transient changes of the real (Aer) (red) and imaginary part (AeJ) (black) of the dielectric 
function of a single 125 nm gold nanoparticle. Pump (probe) fluence is 0.65 J/m^ (0.05 J/m^). 
Dashed lines are calculations as described in the text. 

we can reproduce the position of the absorption cross section (584 nm) obtained including 
retardation using Mie's theory by reducing the plasma frequency to 7.965 eV. However 
when modeling the FWM dynamics (see dashed lines in Fig. HI the particle-dependent pa- 
rameters used in these simulations are Gop = 2.5 x 10^^ W/K and g = 1.8) we can reproduce 
the overall dynamics only qualitatively. We do correctly simulate the signs of AeR_i and 
their relative amplitude, with the positive sign of P being a consequence of probing on the 
blue side of the SPR. However, the delayed rise of Aei observed in the experiment is not 
reproduced by the model. Moreover the proportionality factor ^ between A (P) and Aej 
(Aen) does not scale with the NP volume. We attribute these discrepancies to the failure 
of the dipole approximation. 

V. VALIDITY OF THE ENERGY TRANSFER IN THE TWO-TEMPERATURE 
MODEL 

To demonstrate the validity of the two-temperature model we have investigated the depen- 
dence of the energy transfer rate from the electrons to the lattice Awd/Ar on the difference 
between the electron and lattice temperature — Ti. All these quantities can be retrieved 
directly from the FWM signal when we assume that the total amplitude ^'fwm = yA^ + P^ 
of the field can be approximated as 
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FIG. II. Electron energy transfer rate as a function of T^ — Ti for 38 nm NP with excitation intensity 
of 0.65 J/m^ (black line) at 550 nm and for a 20 nm NP with excitation intensity of 0.60 (red solid 
line) and 1.3 J/m^ (red dashed line) at 550 nm. The dotted lines represent a linear dependence 
with slopes which are proportional to the electron-phonon coupling constant Ggp. 



Efwm = a {t^ - T^) +b{T,- To) . (27) 

The electron and lattice temperature can then be determined from the FWM amplitude 
considering the energy conservation within the NP, which holds for delay times much shorter 
than the typical thermalization time with the environment of some 100 ps, 

7Te 



u 



(Tf-To^) +Ci(Ti-To) (28) 

leading to 

T = (: 

a - 67/(2Ci 



^ pFWM - ^ ^A'^' (29) 



and 

The values of a and b can be found by comparing the dynamics of Tg calculated with the 
two-temperature model and the value of the electron temperature obtained by Eq.|29l The 
value of the energy u has been estimated from the incident energy flux and the nanoparticle 
absorption cross section. Once the values of a and b are determined, the electron energy 
transfer rate can be calculated in a finite difference scheme between the i-th and i + 1-th 



11 



38nm 



1.5 ■ 




11 = 1 

11 = 15 



9- 




1.0 



-0.5 



0.0 



0.5 



10 



100 



T(ps) 



FIG. III. Transient change of the phase of Ae from the ratio between its imaginary and real part, 
measured on the 38 nm NP for an excitation and detection wavelength of 550 nm. Red dashed line 
is a simulation assuming the parameter r/ = 1 in the shift of the plasma frequency (Eq. l32p . The 
solid line assumes r] = 15 as indicated. The black dotted line is the pulse intensity autocorrelation 
measured in the experiment. 

point of the measured FWM amplitude dynamics: 



The values of —Au/At as a function of the electron- lattice temperature difference Tg — Ti 
were plotted in Fig.[TTlfor a 38 nm NP with excitation intensity of 0.65 J/m^ at 550 nm and 
for a 20 nm NP with excitation intensity of 0.60 and 1.3 J/m^ at 550 nm. The data refer to a 
delay range of 0.7-5 ps, i.e. in the cooling regime from the electron gas to the lattice. Within 
the experimental errors, the energy transfer rate depends linearly on Te — T\, confirming the 
validity of the two temperature model. As expected, the amplitude of the slope does not 
depend on the excitation intensity and it is proportional to the value of the electron-phonon 
coupling constant Ggp- 

VI. RATIOMETRIC DETECTION OF THE SIGNAL PHASE 

An intrinsic ratiometric way to represent the measured FWM signal is to plot the phase 
(fi of A? = |Ae| exp{iif) from the ratio between its imaginary and real part (hence A and 
P). This is shown in Fig. lIIII for the 38 nm NP at 550 nm. Noticeably, the visibility of the 
coherent phonon oscillations is improved in ip since it is independent of decay dynamics 
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(31) 
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FIG. IV. Transient change of the real and imaginary part of Ae measured on the 38 nm NP for 
an excitation and detection wavelength of 550 nm, with simulations assuming r] = 1 (dotted lines) 
and 7] = 15 (dashed lines) and for Tb = 70 fs or Tb = 500 fs as indicated. 

common to real and imaginary part of Ae, while a frequency shift in the SPR changes (p. 

We point out that our simulations include the effect of the change in the plasma frequency 
from the excess electron density calculated as: 



The presence of an excess electron density equal to the density of holes generated via the 
interband transitions is corresponding to rj = 1. Indeed all simulations shown in the paper 
for the dynamics of Aej and Ae^ are performed using rj = 1. However, in the dynamics of 
the signal phase we notice that the initial dynamics on the time scale of the excess density 
lifetime (i.e. for r < 0.5 ps) is better reproduced with r/ = 15 (red solid line in Fig. lIIII) . A 
possible physical explanation for 77 = 15 is the influence of the hole density and the creation 
of carriers in specific points of the Brillouin zone with effective masses which are different 
from the average effective mass mo used in the definition of the plasma frequency. 

When looking at the simulations of Aei and Ae^ at 550 nm (see Fig. lIVp the effect of rj 




(32) 
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manifests as a small shift in the initial rise which goes in opposite directions for the real and 
imaginary part of Ae, namely increasing rj shifts Aei toward earlier times and AeR toward 
later times. This effect however does not modify significantly the initial rise time, which 
instead is determined by the d-state Auger thermalization time rt, = 70 fs. To elucidate 
this point we show in Fig. lIVI the effect of changing rt, to 500 fs (equal to the free electron 
thermalization time constant Tf) for both rj = 1 and f] = 15. A longer Tb clearly shows as a 
slower rise time dynamics which cannot be compensated by changing rj and does not agree 
with the measurements. 
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